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Abstract 

Most spherical thin shells, enclosing black body radiation satisfy the domi- 
nant energy condition if they have at least ~ 30% of the total mass-energy. 
Containers with less mass energy, able to sustain high pressures, contain 
mostly unstable radiation. If they have negligible mass energy they are un- 
able to sustain the pressures and the radiation is unstable to gravitational 
collapse. Containers with black holes and radiation in thermal equilibrium, 
considered in the literature, are often unrealistic. 
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1. Introduction 

Thermal equilibrium of black holes and radiation may exist in closed cavities 
(for recent reviews on the subject see papers by Landsberg |I| and Page 
[§]). The walls must be perfectly reflecting and able to resist pressures or 
tensions upon demand. Or so it seems for the resistance of the box is rarely 
questioned. The mass of the wall is ignored or explicitly neglected like in 
Sorkin, Wald and Zhang j3| though far reaching conclusions are sometimes 
drawn. 

From an engineering view point, cavities with properties as asked for in 
black hole thermodynamics are incredible. Forces in the cavity must equi- 
librate the pressure due to radiation, attenuated by self-gravity attraction. 
But radiation pressure is of the order of its mass energy density. Therefore, 
how can a cavity have negligible mass energy (compared to that of radiation) 
and at the same time resist pressure forces comparable to those of radiation 
(divided by the radius of the box R)l Surely some 'energy condition' is 
strongly violated. 

To shed some light on this question we have calculated the mass M s of 
cavities with thin walls (thin spherical shells) and surfaces equal to AtcR 2 . 
The force F on half a sphere satisfies the dominant energy condition Q if 

- M s c 2 /2R <F< M s c 2 /2R. (1) 

The dominant energy condition is the only one to restrain pressures as well 
as tensions. One should not allow for unlimited pressure, this is not physical. 

The results of this calculation are summarized here in broad terms. Finer 
details are given in the text. Let rj be the fraction of total mass-energy in the 
cavity (77 < 1). When rj < 30% there are two different types of equilibrium 
configurations which are both unviable but for different reasons. There are 
cavities very close to the Schwarzschild radius with pressures so high that 
they do not satisfy eq. ([!]). And there are cavities with bigger radii which 
satisfy eq.(|l]) but in which the radiation is either unstable and on the verge 
of gravitational collapse. In practice, realistic cavities with stable radiation 
must have more than 30% of the total mass energy. 

For cavities with rj > 30% there is another novelty. The inverse tempera- 
ture (3 as a function of R has no succession of maxima and minima anymore 
as is typical of equilibrium configurations with equations of state of the form 
P oc p [|. (3 has only one minimum. The energy of the radiation has only 
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one maximum beyond which no equilibrium configuration exists. That maxi- 
mum is also associated with a limit of thermodynamic stability The stability 
limit is smaller for higher 77's. Thus, radiation collapses into a black hole || 
at smaller energies than thought previously [|ll|] . These are the main features 



of 'viable ' cavities, those that are stable and do not violate the dominant 
energy condition. 



2. Elements of the Problem 

We consider spherical cavities with perfectly reflecting walls enclosing black- 
body radiation (photons). The spacetime outside is thus Schwarzschild and 
the metric can be written in terms of the total mass energy M or m = GM 
(with c = 1) as 

ds * = fi _ ^ d e - (l - ^ \ 1 dr 2 - r 2 dQ 2 , r>R (2) 

dtt 2 = d6 2 + sin 2 6 d(f) 2 

R is the 'radius' of the cavity. The radiation inside has pressure P, mass 
energy density p and local temperature T which all vary with r and, in 
Planck units, are related as follows: 

P = i* ,-£1- (3) 

The solution of Einstein's equations is well known, having been studied by 
Klein 0, Chandrasekhar 0, Sorkin, et.al. [Q, Landsberg [Qj and Page [§]. 
Solutions are scale free and all quantities can be expressed in units of R. 
The metric outside can be written as 

rfs 2 = 2u dt 2 dr 2 _ 2 2 

l-2/i(r) ' ~ V ; 

With all clocks referred to the proper time at infinity t, we have a junction 
condition on u(R): 

e MR) = 1 - ^, (5) 

elsewhere, is(r) is given by 

e 4 ^) p (r) = e 4uiR) p(R) (6) 



2 



where 



7T 



2m 



p(R) = — T (R) and W 1 — T(R) = T c 
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R 



(7) 



Given (3 and R, (^) fixes the values for T(R) and p(R). p(r) determines not 
only v{r) but also p(r): 



H(r) 



mir 



mir 



(8) 



p(r) itself is obtained from Einstein's equations. Setting like in Sorkin et al. 

(9) 



dm 2 
q(r) = —— = Anr p and z = mr, 



dr 



one has 



dp 

dz 
dq 
dz 



q-fi, 



2g[l-|g-4/x] 
1 - 2/i 



(10) 
(11) 



From these two equations one can eliminate z and obtain the Tolman-Oppen- 
heimer-Volkoff equation for q(p): 



dq_ 2q[l 
dp 



iff 



4/i) 



l-2/x)( ? - M ) 



with q = p = 



(12) 



This equation is easily integrated; q(p) has been given in other works, but it 
is useful to have q(p) reproduced here (figure 1). With q(p) one obtain q(r) 
by integrating (|TT|) 



R 



exp 



ib (1 - 2/i) dq 
2q(l - \q - Ap) 



(13) 



in which 

q B = AttR 2 Pb (14) 

and pb is defined by R and /3 as we have seen. Finally if we have q(r) we 
obtain p(r) from eq.(^|). u(r) is thus completely defined and depends on (3 
and on 

m , . 

a = -. (15) 



u{r) is independent of the shell, but /x(-R) = m(R)/R depends on the mass 
energy of the cavity. 



3. The Cavities 

Let a be the mass energy density of the thin sphere enclosing the radiation 
and II the pressure (II > 0) or tension (II < 0) in the shell. The relation 
between a, IT and the metrics on both sides are readily found from second 
fundamental forms and from Gauss-Codazzi identities [§]. Following for in- 
stance Goldwirth and Katz PL one finds that 



AwGRa 



and 



AnGRU = - 
2 



l-2/i- y/1 - 2a, fj, = n(R) 
1 - /i+ \q{jj)~ 



1 — a 



(16) 



(17) 



_V1 - 2a ^/\-2\i 

Other useful expressions may be derived in terms of global quantities. First, 
in terms of the mass M s = 47rR 2 a of the shell for (|T^) 

GM X m s 



R 



= Jl - 2fi - Vl - 2a. 



(18) 



Second, in terms of the force F at the rim of half a sphere: 

F = 2irR n, 

equation ([17]) becomes 



GF = f 



a 



VI ~2a y/l-2/j, 



(19) 



(20) 



The dominant energy condition a > and —o < U < o amounts thus to 
a s > and 

-y</<y, (21) 

or in more physical term, to eq. QT|) if c is not set equal to 1. One straightfor- 
ward consequence of equation flTJD and a s > is that 



fi(R) < a. 



(22) 
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4. Equilibrium Configurations and Mass Energy of the 
Shells 

a s does not represent the mass energy of the shell E s which is given by 

Es = fV=9 ^ shell d 3 * = M s . (23) 

Thus, in units of R, 

§ = a s VT^te. (24) 

It follows from (pf) that 
E 



= Vl - 2a 1 - 2/i - VI - 2aJ < Vl - 2a (l - Vl - 2a) < a (25) 
and therefore E s < M. The fraction of mass energy in the shell is defined by 

>,= §<!• (26) 

For a given 77, equation ( |2"1D relates /i to both a and 77. There is some 
advantage in introducing new variables instead of /1 and a: 



y =^\-2^ z = VT^, (27) 

because of eq.([22"D < z < y < 1. Thus eq . (|24]) for n(a,rj) is now replaced 
by an equation for y(z, rf): 

y-lll + U-l),. (28) 



We know from figure 1 that equilibrium configurations exist only for 

/i < /Wr = 0.246 ~ - or y(ii max ) > -j=- (29) 

This limit 1 / V2 divides the equilibrium configurations into two classes de- 
pending on whether the parabolic curve y(z) in eq.(^) has a minimum below 
or above 1 / \[2. The minimum of y is 



Vmir, 



y/vV-v)- (30) 
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The dividing line at rj = r] corresponds thus to yrj (2 — r] ) = A= or 

rjo = 1 - -= ~ 0.293 ~ 0.3. (31) 
v2 

If 77 < 770, then Dmin < y(rjo) (figure 2). Equilibrium configurations exist for 
all values of /x < \i < Umax but in two separate ranges of z. 
The first range 



21 = < z < —= = z 2 32 

is for small values of z. Indeed both Z\ and z 2 are monotonic functions of 
77 and the interval is always close to the Schwarzschild radius, ranging from 
(0,0) for 1] = to (z\ = 0.172,,22 = 0.414) for 77 = which corresponds to 
Ml = 1.03 < < 1.21 = ^. 

2m — 2m — 2m 

The second range of z is 



1 + y-l - 2,,(2 - ,,) 
= ^(2-,) (33) 

and has a big range with j 3 - < 5^- < oo, the lower limit being between 
= 1.21 for 77 = 770 and = 2 for 77 = 0. 

If 77 > 770, 7/ miri is greater than yi^Hmax) and equilibrium configurations do 
not exist for \i greater than some jl(r)) < fi min ~ \. However, the range for 
z is now continuous (see again figure 2): 

z x = — ^— < z<l. (34) 
2-77 ~ ~ v ' 

With 77 > 770, Z\ > 0.172 and ( |3"1| ) corresponds to a total range 1-03 < ^ < 
00, almost down to the Schwarzschild limit. 



5. Linear Series for Radiative Equilibrium 

The three classes of solutions correspond to quite different lines of equi- 
librium configurations. These lines can be deduces from eq. (|7|) and (|) which 
give 

B = (i^) W* = (35) 
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with and z{y,rf) denned by eg. (p8|) and with y = y/1 — 2/i, (|35| ) gives 
rj) at fixed R. 

When T] < T]o, consider first the class of equilibrium configurations for 
which z varies in the interval (z\, z 2 ). Here /i grows from at z = z\ to p ma x 
at z = Z2 and q varies accordingly, as can be seen on figure 1, from point 
O to P. But figure 1 shows also that equilibrium configurations exists for 
decreasing z's, decreasing /i's and further decreasing q([i). Thus for r/ < r] 
and Zi < z < Z2, z and y may oscillate near z 2 a number of times for which 
B{jj) varies just as shown in figure 3a, a result of <?(/i)'s behavior. The class 
of equilibrium for z% < z < 1 are similar to those of figure 3a but in a different 
range of temperature as can be seen in figure 3b where rj = is also drawn. 
The limit 77 — ^ 0, taken by Sorkin et al., amounts to neglect the cavity. Both 
curves of figure 3 are inward spiralling characteristic of equations of state of 
the form P/p =const. |J. 

When rj > 7/0, equilibrium configurations exist only for < p < p(r]) < 
fi m ax- However, when \i varies from to fi, z varies from z\ to ^fzl where 
y = Umin and with increasing value of z, \i decreases again from \x{j]) to zero. 

The corresponding linear series B(fi) is an unwound spiral shown in figure 
4. Here B(fi) tends to infinity twice, for z = Z\ and z — 1, since in both 
limits y — > 1, \i — > and q — > |/x so that q — > and B — > 00. 

6. Viable Cavities 

Among all equilibrium configurations of given 1], what are those that do not 
violate conditions (pT|) ? Or equivalently, in what region of the (y, z) plane or 
the plane or the (fh(R),m) plane does (|21| ) hold for < rj < 1? If 

in equation (|2l|) we replace / as given by eq. (|20D and a s by eq. (|I8|) and use 
(y, z) coordinates defined in ( |27| ) we obtain the following inequalities: 

These inequalities tell us also that 

(i - 5y) (y - z) < ^ < Q + 3j/) (y - z). (37) 

The left hand side is a limit on pressures (n < a), the right hand side on 
tensions (n > — a). 
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6.1. Limits on the pressures 

Normal cavities must resist the radiation pressure. Pressure turns into ten- 
sions when self-attraction of the shell and radiation becomes important. 
Therefore, ordinary cavities are submitted to tension while pressure in the 
wall of the cavity only appear in extremely relativistic case. 

The left hand inequality (|3"7|) corresponding to II < a can be written with 
y replaced by (|28|) and reads then 

2 ~ v) ~ (2 - VV] (l - z 2 ) < (38) 



z 2 > (39) 



Since q > 0, fl3"E| ) always holds if the left hand side is negative, that is, 

2 

5_ 

2-rj 

(|39"D is certainly satisfied for 77 > | = 0.4, but for a given 77, z > Z\ — T]/(2—r]). 
As a result, ( j3~9"D will hold for all z if it holds for z±, that is, when 77 > | which 
is of the same size as ^0 — 0.3. Thus for r] Q > | the pressure in the wall of 
the cavity is never too high. 

Another piece of information follows from (|39D for y — > 1 , that is, [l — > 0, 
and z -> 1 or to 21 (see figure 2). Since for \i <C 1, q — 3/i, one also has 

<?- 3/^ = ^(1-2/)- 3(1-7/), (40) 

because 1 + y ~ 2. We now replace g by ( f40"D into (|38|): 

1. For « -»• 1, M 0, @ becomes 

- 277(1 - z )<(l-j / ). (41) 

This is always satisfied. The pressure is thus never too high in big 
cavities (R > M). 

2. For z —>■ Z\ =77/(2 — 77), R/2m is small and (|38|) becomes 

i- 77 )(l-7 7 )<^(2-7 7 )(l- 2/ ) (42) 

which is never satisfied when 77 < | and y — > 1. The pressure sustained 
by the cavity becomes thus always too high for R —>■ Ri when 77 < 1/3. 
( |4"2"| ) is always correct for 77 > |, but this we already know. 
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6.2. Limits on tensions 



Limits on tensions are obtained from the right hand side of fl3Tp. With y 
replaced by (|28|) one has 



8q < 1 



~ + n ) + (2 - t)) z 2 



l-z 2 ). 



9r] z 

When y — > 1 (// — > 0) and g ~ 3(1 - (see (gD|) ) , then 



(43) 



1. For z — > 1, M — > 0, (43) holds only if 77 > ~. Thus when r; > | ~ r] , 
tensions are never too high. However, for 77 < | and thus also for shells 
with negligible mass energy, tensions are always too high in big enough 
cavities. This means that even when gravity is weak there is always a 
non negligible lower limit to the mass energy of the cavity. 



2. For z 



3f) becomes 



z\) < 



2(1 + 77) 

(2-77)2- 



(44) 



Since z — z\ — ► 0, this inequality is always satisfied and tensions are 
never too big as R — ► R±. 



The limits on tensions or pressures correspond thus approximately to the 
dividing line between the 2 types of solutions (77 < r] or 77 > r/ ). The 
class of cavities with 30% or more mass energy in the shell can resist any 
pressure or tension in all admissible equilibrium configurations. On the other 
hand tensions are always too high in light cavities (rj < 30%) with negligible 
self-gravity (z — > 1). Cavities with negligible mass-energy explode. And 
pressures are always too high in light cavities close to the Schwarzschild 
limit. Such cavities collapse. Figure 5 shows the region in the (a,/i) plane 
in which [ IX j < a (see also figure 3). It also shows the limited domain of 
existence of the 2 classes of equilibrium for rj < i] . It is clear that most of 
these equilibrium configurations have no viable cavities. 
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7. Stability Limits 

Thermodynamic stability limits for spherical perturbations are reached when 
the radiation energy E r has a turning point in a linear series of equilibrium 
configurations. If the mass energy and field energy outside of the cavity are 
negligible, E r /R ~ ~ a, configurations are stable for < ji < Umax = 
0.246 (see figure 3b)[|TTJ. /i max is also a limit where dynamical instability 
for radial perturbations sets in O. All the other configurations of a coun- 
terclockwise spiralling (3{E r ) curve are unstable [jTCfl . A detailed discussion 
of thermodynamic stability, fluctuations and phase transitions in relativistic 
radiation with rj = has been given in Parentani et al. [DJ. Here we are 
interested to know how stability limits change when one takes account of the 
mass in the cavity (rj ^ 0) and the field energy outside. 

Speaking of field energy, we are aware that the subject is controversial 
and, of course, we take sides. E r is surely not M — M s because M — M s > 
even when there is no radiation. The difference must be in the field energy 
outside the cavity Ef (r > R). To calculate Ef, we use Lynden-Bell and 
Katz's method [12]. We replace the cavity with the radiation by a new 



empty cavity of the same radius R and enough mass energy E' s {^ E s ) to 
produce the same gravitational field at r > R. Since the total mass energy 
is M and the mass energy of the shell is E' s , the field energy Ef — M — E' s 
because the flat interior of the new cavity has no energy. Ef/R depends on 
Mj R or a only 



Ef_ 
R 



Thus, radiation energy or 




E r 
~R 



-(l-z) 



E s E 

a 

R R 



f 



(45) 



(46) 



with E s /R given in (p5|) and y, z defined in (|27|), one can write 



e r = z(l-y) = -(2- V )(l-z)(z- Zl ). 
e r is zero at z — Z\ and z = 1 and has one maximum e rmax 



at 



z' = — — — - where e r 



(I-??) 2 

2(2-77)- 



(47) 



(48) 
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The (3{e r ) curves are inward spirals for 77 < i]q or U-type like curves for 
7] > T]o similar to /3(/a) shown in figure 3 and 4. Assuming radiation is 
thermodynamically stable at low temperature (j3 — > 00) all configurations 



for which < e r < e rmax and j3 has its lowest value are stable [10]. The spiral 



of (3{e r ) at higher values of (3 represent unstable configurations. The limits 
of stability in the plane depend on the value of 77. 

We consider first rj > rj for which Z\ < z < 1 and the y(z) line for 
7] > 770 in figure 2. Starting from z = 1 (/i = 0, zero mass energy), we 
pump in slowly energy, evolving through a succession of quasi-equilibrium 
configurations. This decreases both z and y. As we know, y reaches a 
minimum (or 11 a maximum) at z = yfz\ < 1. But yfzi < — jp* = z'\ the 
radiation energy reaches thus its maximum for \i = fi' < jl < \i max . At that 
point obtained from y(z') = y', 

\i = \^-vf-\{l-rif (49) 

and since 77 > 770, 

,'<i(l-^ -1(1-^ = 1-1 = 1 (50) 

Thus, in a linear series of increasing energy, the radiation becomes unstable 
at fi = fjf before reaching fi = Jx{j]) < fi ma x, the higher the mass energy of 
the cavity, the smaller fi'. The limit of stability in the /3(/i) plane is shown 
in figure 4. 

Next consider rj < rjo with z% < z < 1 for which there is a small region 
in the (a,/x) plane where the dominant energy condition holds (see figure 5). 
When 0.235 < 77 < r] Q , z' is greater than Z3 and e r reaches its maximum for 
a fi' < fimax- For < 77 < 0.235, z 3 is greater than z' and e r reaches its 
maximum for fi' = fi max . The limits of stability, where e r is maximum, are 
shown in figure 3a. 

Finally, in the case where 77 < 770 but z\ < z < z 2 the maximum of e r is 
at z = z 2 where /i = fi ma x because Z2 is always smaller than z'. The limit of 
stability is also shown in figure 3b. 



8. Concluding Remarks 

It is clear that most cavities containing thermal radiation and having less 
than ~ 30% of the total mass energy are not realistic in the sense that they 
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have to sustain too high tensions or pressures. Those cavities, able to sustain 
such tensions or pressures contain, however, in most cases, radiation that is 
thermodynamically unstable and would collapse to form a black hole|]]. In 
particular, all cavities with negligible mass energy (rj = 0) that are viable 
contain unfortunately unstable radiation as can be seen in figure 3b. Thus, 
the point raised by Sorkin et al.0 about total entropy becoming infinite 
in cavities with increasing radius or even reaching the Bekenstein [jn| limit 
S = 2ttRM does not apply to stable radiation in 'viable' cavities. 

One can have some fun by putting cavities with radiation on the other 
side of the Einstein-Rosen bridge in a 'perpetual' Schwarzschild spacetime. 
In this case, a is higher and II is always negative and instead of (|16D and 
( |TTD one has now 



AnGRa = yjl - 2fi + y/l-2a (51) 

and 

l I- a | l-// + g(/x)/3 j 



AuGRTl = — 
2 



Viable cavities are for II > —a. One can also put a small black hole in 
the center like Hawking |TJJ] and look for the changes our considerations 



would make in a thermodynamic stability analysis (Parentani et al.[ll|) if one 



wishes to be realistic about probabilities of phase transition, i.e. probabilities 
for black body radiation to collapse and form a black hole 0. 
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Fig.l: This represents the solution q(fx) of the TOV equation (12) and 
/2 max ~ 0.246. 

Fig. 2: Two y(z) lines, parametrized in /i: one for rj = 0.35 > r/ ~ 0.29 
whose minimum y min is at z = \J zi(rj) for fj, — ft, < fi max and another for 
7] = 0.15 < r] ~ 0.293 whose minimum at z = \J zi(rj) is below y m i n (vo)- 

Fig. 3a: The dark region in this -B(/i) plane corresponds to |II| < a when 
z\ < z < z 2 and holds only for 0.25 < i] < r] ~ 0.29. There are no viable 
cavities for rj < 0.25. The two linear series in this figure are in solid lines 
where configurations are stable (section 7) and dashed for unstable ones. 

Fig. 3b: Dark regions in the plane correspond to |n| < a when z% < z < 1 
and holds for < rj < i] . the two linear series in the figure are in solid lines 
where stable (section 7) and dashed where unstable. There are no stable 
lines in the dark region for r/ < 0.11. The dot-dashed line represents i] = 
configurations. As can be seen, such cavities either contain unstable radiation 
or are unable to sustain the pressures. 

Fig. 4: The linear series B(fi) for rj = 0.35 is characteristic of all lines 
for rj > i] . The dashed part of the line represents unstable configurations 
(section 7) the thick lines are for stable ones. All configurations are unstable 
above the thin continuous line representing the limit of stability. 

Fig. 5: The dark region of the (a,/i) plane are points associated with 
equilibrium radiation in cavities satisfying the dominant energy conditions. 
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